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ABSTRACT 

Spectral and kinematic studies suggest that the nonthermal radio source Sgr 
A*, located at the center of the Milky Way, is a supermassive compact object 
with a mass ~ 2 — 3 x 1O 6 M . Winds from nearby stars, located ps 0.06 pc to 
the east of Sgr A*, should, in the absence of any outflow from the putative black 
hole itself, be accreting onto this object. We report the results of the first 3D 
Bondi-Hoyle hydrodynamical numerical simulations of this process under the 
assumption that the Galactic center wind is generated by several different point 
sources (here assumed to be 10 pseudo-randomly placed stars). Our results 
show that the accretion rate onto the central object can be higher than in the 
case of a uniform flow since wind-wind shocks dissipate some of the bulk kinetic 
energy and lead to a higher capture rate for the gas. However, even for this 
highly non-uniform medium, most of the accreting gas carries with it a relatively 
low level of specific angular momentum, though large transient fluctuations can 
occur. Additionally, the post-bow-shock focusing of the gas can be substantially 
different than that for a uniform flow, but it depends strongly on the stellar 
spatial distribution. We discuss how this affects the morphology of the gas in 
the inner 0.15 pc of the Galaxy and the consequences for accretion disk models 
of Sgr A*. 

Subject headings: black hole physics — hydrodynamics — Galaxy: center- 
galaxies: nuclei — ISM: jets and outflows — stars: mass-loss 
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1. Introduction 

Sgr A* may be a massive (~ 2 — 3 x 10 6 M ) point-like object dominating the 
gravitational potential in the inner ^ 0.5 pc region of the Galaxy. This inference is based on 
the large proper motion of nearby stars (|Haller, et al. 1995| ; [Eckart fc Genzel 19971 Oenzel 



Eckart, Ott fc Eisenhauer 1997|) , the spectrum of Sgr A* ( [Melia, Jokipii, fc Narayanan 



1992|) , its low proper motion ( ^ 20kms 1 ; [Backer 1996| ), and its unique location ( [Lacy. 



[Achtermann, fc Serabyn 1991[ ). The gaseous flows in this region are themselves rather 



complex, and key constituents appear to be the cluster of mass-losing, blue, luminous stars 
comprising the IRS 16 assemblage, which is located within several arc seconds (1" ~ 0.04 
pc in projection at the distance to the Galactic center) from Sgr A*. Measurements of 
high outflow velocities associated with IR sources in Sgr A West ([Krabbe, et al. 19911 ) and 
in IRS 16 QGeballe, et al. 1991|) , the H2 emission in the circumnuclear disk (CND) from 



molecular gas being shocked by a nuclear mass outflow flGatley, I., et al. 1986|) , broad Bra 
Br7 and He I emission lines ( [Hall, Kleinmann fc Scoville 1982| ; |Allen, Hyland fc Hillier 1990 



Gcballe, et al. 1991|) , and radio continuum observations of IRS 7 ( [Yusef-Zadch fc Melia 



|1991| ), provide clear evidence of a hypersonic wind, with velocity v w ~ 500 — 1000 kms -1 , a 
number density n w ~ 10 3-4 cm -3 , and a total mass loss rate M w ~ 3 — 4 x 10~ 3 M Q yr -1 , 
pervading the inner parsec of the Galaxy. Many of Sgr A*'s radiative characteristics may 
be due to its accretion of the IRS 16 wind. 



In the classical Bondi-Hoyle (BH) scenario ( [Bondi fc Hoyle 1944! ), the mass accretion 



rate for a uniform hypersonic flow is Mbh = KRA 2 mHn w v w , in terms of the accretion radius 
Ra = 2GM/v w 2 . At the Galactic center, with n w ~ 5.5 x 10 3 cm -3 and v w ~ 700 kms -1 , 
we would therefore expect an accretion rate Mbh ~ 10 22 gs _1 onto the black hole, with a 
capture radius Ra ~ .02 pc. Since this accretion rate is sub-Eddington for a one million 
solar mass object, the accreting gas is mostly unimpeded by the escaping radiation field 
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and is thus essentially in hydro dynamic free-fall starting at Ra- Our initial numerical 
simulations of this process, assuming a highly simplistic uniform flow past a point mass 
( |Ruffcrt fc Melia 1994j ; |Coker fc Melia 1996|) have verified these expectations. 



On the other hand, the nature of Bondi-Hoyle accretion onto a point like object also 
presents somewhat of a challenge in understanding what happens to the gas as it settles 
down into a planar configuration close to the event horizon. Fluctuations beyond the bow 
shock (located at ~ Ra) produce a transient accretion of net angular momentum that ought 
to result in the formation of a temporary (albeit small) disk. The circularization radius is 
r c « 2X 2 r g , where r g = 2GM/c 2 is the Schwarzschild radius and A is the specific angular 
momentum in units of cr g . Our simulations of the BH accretion from a uniform flow suggest 
that in this configuration (A) ~ 3 — 20. More realistically, the inflow itself carries angular 
momentum, so that the formation of a disk-like structure at small radii (i.e., r rs 10 2 ~ 3 r g ) 
may be difficult to avoid ( Melia 1994|) . However, the observations do not appear to favor 



the presence of a standard a-disk at small radii ( Melia 1996| ). The current upper limits on 



the infrared flux from Sgr A* (Menten, et al. 1997|) suggest that either (1) the circularized 



flow does not form an a-disk, but rather advects most of its dissipated energy through the 
event horizon (e.g., Narayan, Yi, fc Mahadevan 1995| ), (2) the Bondi-Hoyle flow merges 



into a massive, fossilized disk, storing most of the deposited matter at large radii ( |Falcke & 
Melia 1997|) , or (3) Sgr A* is not a point-like object ( [Haller, et al. 1995| ). 



Added to this is the fact that in reality the flow past Sgr A* is not likely to be uniform. 
For example, one might expect many shocks to form as a result of wind-wind collisions 
within the IRS 16 comples, even before the plasma reaches Ra- With this consequent loss of 
bulk kinetic energy, it would not be surprising to see the black hole accrete at an even larger 
rate than in the uniform case. The implications for the spectral characteristics of Sgr A*, 
and thus its nature, are significant. We have therefore undertaken the task of simulating 
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the BH accretion from the spherical winds of a distribution of 10 individual point sources 
located at an average distance of a few Ra from the central object. As we shall see below, 
the accretion rate depends not only on the distance of the mass-losing star cluster from the 
accretor but also on the relative spatial distribution of the sources. As suggested by related 
work involving linear gradients ( [Ruffert fc Anzer 1995| ), the average value of A is also larger 
than that of the uniform case, exhibiting large temporal fluctuations, but still not as large 
as one might expect. In this Letter we present the results of the first 3D hydrodynamical 
calculations of multiple- source stellar winds accreted by a point mass. 



2. The Calculational Algorithm and Physical Setup 

Although we here attempt to model the wind source more realistically than in previous 
work, we shall still focus on the hydrodynamical aspects of the flow and thus we do not 
include heating or cooling (see, e.g., Melia 19941 ). Nor do we include any relativistic effects. 



Instead, we assume the medium to be an unmagnetized, polytropic gas with pressure 
P = (7 — l)pe where p, and e are the mass and internal energy densities, respectively, 
and 7 = 4/3 is the adiabatic index. This value of 7 was chosen so that the sonic point 
d s = Ra x (5 — 37) /4 occurs outside the accretor, which typically is modeled with a radius 
~ O.IRa- Some earlier work (e.g., Petrich, et al. 1989]) found that varying 7 affects the flow 



pattern significantly, particularly at radii much smaller than Ra, while others (e.g., [Ruffert 
1996) found that it does not, particularly at the larger radii considered here. 



We use the numerical algorithm ZEUS-3D, a general purpose code for MHD fluids 
developed at NCSA ( Norman 1994 ). The calculations are carried out within a cubical 



domain of solution with 112 3 active zones geometrically scaled so that the central zones 
are 1/32 times the size of the outermost zones, closely mimicking the "multiply nested 
grids" arrangement used by other researchers (e.g., |Rutfert fc Melia 1994j ). This allows 
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for maximal resolution of the accretor within the computer memory limits available while 
sufficiently resolving the wind sources and minimizing zone-to-zone boundary effects. The 
total volume is (16Ra) 3 or ~ (0.28 pc) 3 using a point mass M — 1 x 10 6 M Q located at the 
origin and a Mach 10 Galactic center wind with v w = 700 kms -1 . 

The initial density is set to a small value and the velocity is set to zero. The internal 
energy density is chosen such that the temperature is ~ 10 4 K. Free outflow conditions 
are imposed on the outermost zones and each time step is determined by the Courant 
condition with a Courant number of 0.5. The 10 identical stellar wind sources are modeled 
by forcing the velocity in 10 subregions of 5 3 zones to be constant (at v w ) while the densities 
in these subvolumes are set so that the total mass flow into the volume of solution, M w , is 
3xl0~ 3 M yr _1 . Recent observations flNajarro, et al. 1994[ ) suggest that the IRS 16 wind 



is somewhat colder than the temperature assumed here, with a Mach number of ~ 30, and 
that it may originate from more than 15 sources of varying strength rather than 10 sources 
of equal strength. This level of detail will be addressed in future work. 

To gauge the dependence of our results on the source configuration, we here simulate 
the BH process using two different stellar distributions. For numerical reasons, the sources 
need to be placed a minimum of 8 zones ( ^ -2Ra) apart so no attempt is made here to 
imitate exactly the distribution within the IRS 16 cluster. The average location for the 
sources in run 1 is ~ AR^z (with the accretor located at the origin) and ~ 3Raz for run 2. 
For run 1, the sources are distributed fairly randomly with x, y, and z being allowed to vary 
by ±3Ra of the average location. For run 2, x and y are allowed to be as much as 5 Ra 
from the average location while z can only vary by ±1Ra- These distributions are chosen to 
represent the extremes of a spherical stellar distribution versus a more or less planar one. 
For comparison, IRS 16 N W and IRS 16 NE, two bright members of the IRS 16 cluster, are 
~ 3 and ~ 7Ra away (in projection) from Sgr A* ( Menten, et al. 1997] ) . Each simulation is 
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evolved for 2000 years. The wind crossing time is ~ 400 years and equilibrium, the point at 
which the original gas is swept out of the volume of solution and the mass accretion rate 
stabilizes, appears to be reached within 2 crossing times, or ~ 800 years. 

3. Comparison of the Flow Patterns 

Figure 1 shows a logarithmic negative grey scale image of the density profile for a slice 
running through the center of the accretor, for run 1 taken 2000 years after the winds are 
"turned on". The image is 6 Ra on a side with the O.IRa radius accretor located at the 
center. One stellar source can be clearly seen to the upper right of the accretor; the other 
9 sources are either not in the plane of the slice or are off to the right. The image spans 4 
magnitudes of density with white being £ 0.01n w and black being £ 10n w . Figure 2 is a 
similar image for run 2. Note the large density fluctuations across the overall region. Also, 
the flow pattern is clearly different for the two configurations. However, once the stellar 
winds have cleared the region of the original low density gas, both simulations point to an 
overall average density (~ 10 3 cm -3 ) in agreement with observations. 

A substantial difference between these point source simulations and previous uniform 
calculations is that although there are transient fragmentary shocks extending from the 
accretor boundary (see Figure 1), there is no large-scale bow shock structure. There 
are, however, bow shocks around some of the stellar sources as their wind impacts the 
stronger accumulated wind of the other sources (see the source in Figure 1), and, for run 
2, a semi-stable partial bow shock appears in front of the accretor (see Figure 2). The 
wind-wind collisions convert kinetic energy into thermal energy with the result that a larger 
fraction of the gas is captured by the central engine. Although this effect is likely to be 
sensitive to the actual stellar distribution, a larger density of point sources should produce 
more wind-wind collisions and perhaps further raise the accretion rate. 
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4. The Accretion of Mass and Angular Momentum onto the Central Object 

For a wind originating from a single point source located at a distance R (> Ra) from 
the accretor, M is less than M BH due to the divergent flow. Specifically, for R = AR A) 
M ~ 0.15 M BH . However, wind-wind collisions from multiple sources reduce the effective 
v W) thereby increasing the effective Ra so that all of the shocked gas will then accrete and 
more than compensate for the geometric losses. In Figure 3 we present the mass accretion 
rate, M, and the accreted specific angular momentum, A, versus time for run 1, starting 2 
crossing times (~ 800 years) after the winds are "turned on". Figure 4 presents the results 
for run 2. The average values for the mass accretion rate once the systems have reached 
equilibrium are M = 2.1 ± 0.3Mbh for run 1 and 1.1 ± 0.2Mbh for run 2. 

The mass accretion rate shows high frequency temporal fluctuations (with a period 
of ^ 0.25 yr) due to the finite numerical resolution of the simulations. The low frequency 
aperiodic variations (on the order of 20% in amplitude) reflect the time dependent nature 
of the flow. Thus, the mass accretion rate onto the central object, and consequently the 
emission arising from within the accretor boundary, is expected to vary by ^ 20 — 40% 
(since in some models L may vary by as much as oc M 2 ) over the corresponding time scale 
of < 100 years, even though the mass flux from the stellar sources remains constant. 

Similarly, for run 1, the accreted A can vary by 50% over ^ 200 years with an average 
equilibrium value of 37 ± 10. For run 2, (A) = 62 ± 5 with smaller amplitude long term 
variability. Since previous uniform simulations resulted in (A) ~ 3 — 20, it appears that 
even with a large amount of angular momentum present in the wind, relatively little specific 
angular momentum is accreted. This is understandable since clumps of gas with a high 
specific angular momentum do not penetrate within 1 Ra- The variability in the sign of 
the components of A suggests that if an accretion disk forms at all, it dissolves, and reforms 
(perhaps) with an opposite sense of spin on a time scale of ~ 100 years. 
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5. Discussion, Conclusions, and Future Work 

A variety of accretion scenarios for Sgr A* have been proposed over the years flMclia 



1994| ; |Kalcke fc Heinrich 1994] |Narayan, Yi, fc Mahadevan 1995| ; [Beckert, et al. 1996Q , 



each with its own restrictions on M and A. While the accreted specific angular momentum 
determined in the present simulations is an order of magnitude too small to support the 
fossil disk scenario (since then the energy liberated as the wind impacts the fossil disk 
should be visible; |Falcke fc Melia 1997| ), it is still large enough that any standard a-disk 
would be easily detectable ( |Melia 1994| ). The advection dominated disk scenario ( [IN arayan 



Yi, fc Mahadevan 1995), while permitting a large range of values for A, requires an accretion 



rate of ^ 10~ 5 M^yr" 1 or roughly O.IMbh, which is 10-20 times smaller than the value 
derived here. In addition, our M is more than 10 times larger than that permitted by 
the "mono-energetic" electron model of Sgr A* QBeckert, et al. 1996|) . It appears that 



additional work is needed to reconcile disk models with the fact that the observed multiple 
wind sources result in a large mass accretion rate onto the central engine, if its mass is 
~ 10 6 M Q . In view of this, it may not be unreasonable to conjecture that in fact no flattened 
disk actually forms in Sgr A*, but rather that the excess angular momentum is dissipated 
in a quasi-spherical configuration and that the thermalized energy is then advected inwards 
through the event horizon before the gas settles onto a plane ( |Melia 1992| ). 

These simulations suggest that the ~ 0.1 pc region of the Galaxy, centered on the wind 
sources, is swept clear of gas, leaving a hot, low density gas filling the central cavity. This 
is consistent with observations of the region within the CND ( |Yusef-Zadeh fc Wardle 1993| ) , 
and may be acting in concert with other mechanisms to produce the sharp inner edge of 
the CND (e.g., the abrupt change in gravitational potential; |Duschl 1989| ). Additionally, a 
tongue of hot, dense gas has been observed that connects members of the IRS 16 cluster to 
Sgr A* QGeballe, et al. 199 1| ) . It is worthwhile noting that the images in Figures 1 and 2 
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show ridges of dense gas connecting the sources in the figures to the accretor. 
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Fig. 1. — A logarithmic negative grey scale snapshot, 6 Ra on a side, of the gas density for 
run 1 taken 2000 years after the winds are "turned on" . It is a slice in the x — z plane (with 
z to the right) through the accretor (marked with a filled white circle). White corresponds 
to a number density & 0.01n w and black corresponds to a number density £ 10n w . 

Fig. 2. — Same as Figure 1 but for run 2 at 2000 years. 

Fig. 3. — Accretion results for run 1. The upper solid curve is accreted specific angular 
momentum A (in units of cr g ). The scale for A is on the left side. The lower dotted curve 
is the mass accretion rate M (10~ 4 M Q yr _1 ) versus time. The scale for M is shown on the 
right side. 

Fig. 4. — Accretion results for run 2. 
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